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ABSTRACT 

We present numerical time-dependent calculations for fall-back disks relevant to GRBs in which the disk of 
material surrounding the black hole (BH) powering the GRB jet modulates the mass flow, and hence the strength 
of the jet. Given the initial existence of a small mass < 10~^Mq near the progenitor with a circularization 
radius ~ 10^'^ — 10^^ cm, an unavoidable consequence will be the formation of an "external disk" whose outer 
edge continually moves to larger radii due to angular momentum transport and lack of a confining torque. For 
long GRBs, if the mass distribution in the initial fall-back disk traces the progenitor envelope, then a radius 
^ 10^^ cm gives a time scale ^ 10'* s for the X-ray plateau. For late times t > 10^ s a steepening due to a 
cooling front in the disk may have observational support in GRB 060729. For short GRBs, one expects most 
of the mass initially to lie at small radii < 10^ cm; however the presence of even a trace amount ~ 10~^Mq 
of high angular material can give a brief plateau in the light curve. By studying the plateaus in the X-ray decay 
of GRBs, which can last up to ~ 10^ s after the prompt emission, Dainotti et al. find an apparent inverse 
relation between the X-ray luminosity at the end of the plateau and the duration of the plateau. We show that 
this relation may simply represent the fact that one is biased against detecting faint plateaus, and therefore 
preferentially sampling the more energetic GRBs. If, however, there were a standard reservoir in fall-back 
mass, our model can reproduce the inverse X-ray luminosity-duration relation. We emphasize that we do not 
address the very steep, initial decays immediately following the prompt emission, which have been modeled 
by Lindner et al. as fall-back of the progenitor core, and may entail the accretion of > IMq. 

Subject headings: Accretion, accretion disks — Gamma ray burst: general — Gamma-ray burst: individual: 
GRB 060729, GRB 05 1221 A 



1. INTRODUCTION 

The tentative pie-Swift hints that breaks occur in the long- 
term afterglow light curves simultaneously in different wave- 
bands at - 1 - 10 d (Frail et al. 2001) after the GRB have 
not been borne out by numerous detailed observations of Swift 
GRB afterglows in recent years (Ghisellini et al. 2007, Gates 
et al. 2007, Racusin et al. 2008, Liang et al. 2008). With 
Swift (Gehrels et al. 2004), the long-term (t > 10^ s) X-ray 
behavior has turned out to be surprisingly complex, exhibit- 
ing alternating steep and shallow slopes when plotted as log 
Fx versus log(t — TO), where TO is the trigger time for prompt 
BAT emission (Zhang et al. 2006; Nousek et al. 2006). 

The traditional viewpoint that long-term fading of the GRB 
afterglow represents the deceleration and spreading of a rela- 
tivistic jet (e.g.. Sari, Piran, Halpern 1999, Frail et al. 2001) 
has given way to a much broader effort involving a diverse 
variety of hypotheses by which a central engine can be pow- 
ered for a long time, including such scenarios as magnetars 
(Fan & Xu 2006; Toma et al. 2007; RowHnson et al. 2010), 
quark stars (Staff, Niebergal, & Guyed 2008), and long term 
accretion onto a BH formed during the collapsar or NS-NS 
merger (Kumar et al. 2008a, 2008b; Metzger et al. 2008a, 
2008b). The accretion hypothesis entails exploring the pos- 
sibility that the long-term decay of the X-ray flux is not due 
to the deceleration of baryonic ejecta, but rather a secular de- 
crease in the rate of accretion powering the central engine, 
and therefore indirectly the jet (Kumar et al. 2008a, 2008b; 
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Metzger et al. 2008a, 2008b, Cannizzo & Gehrels 2009, here- 
after CG09). For long GRBs (IGRBs), the early, steep rate of 
decay may be giving us information about the radial density 
distribution within the progenitor core (Kumar et al. 2008b), 
whereas the later decay may be governed by the outward ex- 
pansion of the transient disk formed from the remnants of the 
progenitor (CG09). As regards short GRBs (sGRBs), even if 
a small amount of material (~ 10~^ — 10"'' M©) is expelled 
during the NS-NS merger and later accreted in a disk, that 
would be sufficient to power a bright afterglow, which may 
also be strongly influenced by the effects of r-process nucle- 
osynthetic heating in the neutron rich material that becomes 
the disk (Metzger et al. 2010). 

Zhang et al. (2006) present a schematic for the decaying 
GRB light curve as seen by the XRT on Swift. The decay 
is traditionally shown in log F — log t. There are four basic 
power-law decay {F oc i~") regimes: (i) a steep decline fol- 
lowing the prompt emission with ai ~ 3 out to 10^ — 10"^ s, 
(ii) a plateau with an ~ 0.5 out to 10^ — 10^ s, (iii) a steep- 
ening with am ~ 1.2 out to 10^ — 10^ s, and (iv) a further 
steepening at late times (not always seen) with ajv — 2. 

CG09 present a general analytical formalism to explain the 
different power law decays using a fall-back disk, where the 
variations in a could potentially be explained by different 
physics operating within the disk. The results of CG09 were 
purely analytical; in this work we present time dependent nu- 
merical calculations in order to examine in more detail the 
potential of the model, and we apply the results to XRT data 
for one IGRB and one sGRB, taking the best studied of each 
class. In Section 2 we review the Dainotti relation, an empiri- 
cal relation involving the duration and luminosity of segment 
11, in Section 3 we present our detailed numerical model, in 
Section 4 we compare the model with observations for the 
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IGRB 060729, the GRB with the longest observational time 
series in X-rays, in Section 5 we compare theory and data 
for the sGRB 051221A, the sGRB with the longest and most 
detailed XRT light curve, in Section 6 we revisit the Dainotti 
relation, in the context of our numerical results, and in Section 
7 we discuss and summarize our results. 

2. THE DAINOTTI RELATION 

Dainotti et al. (2008, 2010) found an empirical relation be- 
tween the duration of the X-ray plateau in the source frame 
tjj = tu{l + z)~^, and the X-ray luminosity at the end 
of the plateau, also corrected into the source frame, ijj. 
Expressing their relation in the form logj^g Ljj(erg s^^) — 
a + 51ogiQtJj(s), Dainotti etal. (2010)finda = 51.1±1 and 
b = —1.1 ± 0.3. Their sample of IGRBs with small errors is 
defined by w < 4, where u = [(logJLJj)^ + (log (5tjj)^]^/^, 
and log JLjj and log 6t^^ represent the logarithmic errors in 
Ljj and tfi, respectively. Using the same sample of 62 IGRBs 
with small errors taken from Dainotti et al. (2010), we find^ 
a = 50.5 ± 0.4 and b = -0.92 ± 0.1, close to their results. 
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if long-term accretion is the correct explanation for the long- 
term X-ray light curve of GRBs, that would have implica- 
tions for the amount of mass in the initial fall-back disk. To 
first order, the fact that L^^t^^ is constant would imply the ac- 
creted mass reservoir is constant. Willingale et al. (2007) 
present a simple formalism for integrating the X-ray flux on 
the plateau and subsequent am ~ 1.3 decay to obtain a to- 
tal energy for segments II and III. We can convert this into 
a mass by making a few plausible assumptions about the en- 
ergetic efficiencies (Krolik, Hawley, & Hirose 2007), e.g., an 
X-ray afterglow beaming factor / ~ lO^'^ — 10^^, and an effi- 
ciency enct — 0.01 — 0.1 with which accretion onto the inner 
engine powers the observed X-rays, presumed to be created 
within the beamed jet. This overall accretion efficiency has 
two components, (1) the efficiency eacc with which accretion 
onto the BH powers the jet, where Cacc — 0.4 for the antici- 
pated high spin BHs in GRBs, and (2) the efficiency ex with 
which the jet power is converted into 0.3 — 10 keV X-rays 
that can be observed by XRT. Therefore Enot = ^acc^x- We 
use eqns. [3] and [4] from Willingale et al. (2007) to obtain a 
total energy SEx and hence SAI derived from the X-ray flu- 
ence, for segments II and III. The results are shown in Figure 
1 . The hatched area is bounded from above by an estimate of 
the effective limiting XRT flux^ for being able to observe and 
characterize plateaus sufficiently well that they would satisfy 
u < 4. There appears to be a selection effect giving rise to 
the Dainotti relation, namely that the flux detection limit for 
XRT prevents one from clearly observing and parameteriz- 
ing faint plateaus at high z. The Dainotti relation, effectively 
equivalent to ijjijj — constant, results from the sampling of 
GRBs primarily beyond z ~ 1.5, for which faint plateaus 
are observationally biased against. Thus we are seeing basi- 
cally a narrow strip corresponding to the upper end of a much 
broader distribution, which makes it appear that in^n — 
stant. For z < 1.5 one sees curvature in the lower envelope 
of SEx values due to the strong dependence of the detection 
limit on z. Nevertheless, the fact that there appears to be a 
well-defined upper limit SM « 10^'' — 10~'^Mo is interest- 
ing: if one starts with a ~ lOAf© progenitor and if accretion 
governs the long-term X-ray light curve, this indicates that no 
more than a fraction ~ 10~^ — 10"** of the progenitor mass 
survives the hypernova to form a fall-back disk (excluding 
segment I). Given the potential for strong outflows during the 
accretion process, this 6M should be viewed as a lower limit. 



Fig. 1 . — Inferred total accretion mass for the plateau -|- later decay phases 
of GRBs, using data from Dainotti et al. (2010) for their 62 well-constrained 
IGRBs. The hatched area (shown in red) indicates a putative limiting XRT 
detection flux level /n = 10^^^ erg cm^^ s^^ (adopting a plateau duration 
til = 10^ s) for being able to study a plateau to sufficient accuracy that it 
would have been included in the Dainotti et al. sample of 62 GRBs with 
good statistics and known redshifts. (Three of the redshifts given by Dainotti 
et al., z = 0.08 for 051109B, z = 0.78 for 060202, and z = 1.16 for 
0705 IB, are less secure than the others, i.e., based either on photometry or 
X-ray spectroscopy, and have been indicated by blue symbols.) We adopt 
a beaming factor / = 1/300 and a net efficiency for powering the X-ray 
flux enot = eacc^jf = 0.03 to Convert from X-ray fluence (i.e., total X-ray 
energy, after taking into account Airdj^) to accreted mass. 

If the inverse relation between ijj and tjj is physical, and 

' We employ a different method than Dainotti et al. (2010), a Monte Carlo 
technique in which 10® data sets are created with Icr errors randomly either 
added or subtracted to each data point in both x and y directions. Thus for A'^ 
data points one could in principle have Af* distinct data sets. For each data 
set, a and b values are calculated, and the final values for a, 5a, b, and Sb ai'e 
taken from the averages and standard deviations of the lO'' values. 



3. ACCRETION DISK PHYSICS 

By writing the equations for mass continuity and angular 
momentum transport in cylindrical coordinates, assuming Ke- 
plerian rotation — GMbht^^ and integrating over the 
vertical thickness of the accretion disk, one arrives at an equa- 
tion for the evolution of the surface density E = 2ph, where p 
is the density and h the disk semithickness (actuaUy pressure 
scale height). 
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* Our adopted XRT plateau detection level 10~^^ erg cm~^ 
related to the nominal single pointing detection limit ~ lO"^** erg cm~^ 
s~^, and in fact exceeds it considerably: In order to distinguish the differ- 
ent phases of the canonical light curve (i.e., plateau versus steep decay), one 
must take into account the effect of Swift orbital gaps, which places a much 
stronger constraint on the light curve characterization than simple flux detec- 
tion. In addition, the presence of flares may further complicate the picture, 
but they usually are confined to t < 10^ s and generally do not represent a 
major portion of the total energy budget. 
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The kinematic viscosity coefficient 

where P is the pressure and ass is the Shakura-Sunyaev 
parametrization of the angular momentum transport and heat- 
ing (Shakura & Sunyaev 1973). Guided by the current concor- 
dance between inferred values of ass — 0.1 from dwarf nova 
outburst decays, both fast, thermal decays (Smak 1984) and 
slow, viscous decays (Cannizzo et al. 2010), and global 3D 
general relativistic magnetohydrodynamic (GRMHD) models 
of the magnetorotational instabihty (MRI) which also give 
ass — 0.1 at large radii in the disk (McKinney & Narayan 
2007 ab), we set ass = 0.1. Equation (1) is discretized fol- 
lowing the method of Bath & Pringle (1981), in which grid 
points are distributed as r^/^. In addition, one has a thermal 
energy equation governing the temperature evolution, 

dT 2(A-B + C + D) UTl d , , dT 

IJT = TT 

ot Cpi-, fiCp r or or 

where the viscous heating A = (9/8)i^il^E, the radiative 
cooling B = (jTg , and C and D represent radial heat fluxes 
due to turbulent and radiative transport (Cannizzo et al. 2010). 
For the calculations we present in this work, the thermal 
equation only becomes of potential importance at very late 
times > 1 yr. In other words, A — B to very good pre- 
cision for most of the evolution, due to the fact that ( = 
rflogr(E)/dlogE > 0. In other words, disks with C > 
are viscously and thermally stable (Piran 1978), and therefore 
simple viscous evolution provides a satisfactory physical de- 
scription. 

After ~ 1 yr the surface density in the very outer disk drops 
below that associated with the hysteresis due to the transition 
between ionized and neutral gas, the source of the dwarf nova 
limit cycle mechanism. This launches a cooling front that 
propagates to smaller radii and rapidly shuts off the supply 
of gas to the inner disk, as the effective viscosity plummets. 
The dwarf nova hmit cycle is an accretion disk instabihty in 
which the T'(S) relation exhibits a hysteresis at a temperature 
corresponding roughly to the peak in the Rosseland opacity 
curve — effectively an "S"-curve when plotted as T versus 
E. It accounts for the outbursts seen in dwarf novae (Lasota 
2001). The global manifestation of the instability is that ma- 
terial accumulates in the disk during quiescence and accretes 
onto the central star during outburst. The transition between 
the two states is mediated by the action of heating and cooling 
fronts. The cooling front we see at late times in our calcula- 
tions are a result of the disk mass becoming so small that the 
transition from ionized to neutral gas is instigated. (A detailed 
discussion of the effects associated with the heating and cool- 
ing fronts is given in Cannizzo [1993, 1998].) 

At early times and small radii the local disk accretion can be 
significantly super-Eddington, and in the thin disk formalism 
one would have h/r >> 1, an inconsistency. Section 2.3.2 of 
CG09 contains scaUngs for the slim disk branch of solutions 
(Abramowicz et al. 1988) which we employ when h/r > 
(3/4)1/2 locally. Combining eqns. (24) and (26) of CG09 
gives 



T(S) = 1.56 X 10^ K mB/n.ir-n 



other types of disks are possible, such as neutrino-cooled ac- 
cretion disks (e.g., Narayan, Piran, & Kumar 2001; Kohri, 
Narayan, & Piran 2005; Chen & Beloborodov 2007). The 
controlling dynamics of the long term accretion rate is gov- 
erned by the physical state of the outer disk, since that is 
where the slowest (i.e., controlling) time scale is. A com- 
plete physical modeling of the disk at small radii, including 
neutrino coohng and nucleosynthesis, is beyond the scope of 
this work, but also not necessary given in our primary inter- 
est, the long term accretion rate. At late times and large radii 
where h/r < (3/4)^/^, we utilize the thin disk scalings from 
CG09. Combining eqns. (10) and (1 1) from CG09 gives 



T(S) = 84.1 K ag' 



1/6 

'BH,!' 11 



'l/2y.2/3 



(5) 



(4) 



where mBH,i = Mbh/IIOM©) and rn = r/(10" cm). 
For completeness, we note that at early times and small radii. 



where ass,-i = ass/0-1. From a pragmatic standpoint, even 
though the use of a thin disk at small radii would be incon- 
sistent because one would find h/r » 1, the use of a slim 
disk at small radii versus a thin disk has minimal effect on the 
Ught curve. This comes about because the controlUng time 
scale in the problem comes from the largest disk radius. In 
fact, one could probably model the entire disk as a one-zone 
model, taken at the outer disk edge as was done by Metzger 
et al. (2008a) for NS-NS merger disks, and still get a reUable 
result. 

The composition of the accretion disk is substantially dif- 
ferent from solar. The T(E) relation exhibits a weak in- 
verse scaling with mean molecular weight ji for gas pres- 
sure and electron scattering opacity, T cx iJ,~^^^ (Cannizzo & 
Reiff 1992). The expected mixture of alpha elements yields 
/X ~ 16, rather than the value /z ~ 0.65 more generally rele- 
vant for solar composition material. Thus T cx /i^i/'' brings 
about a decrease in T(E) by a factor ~ 2.9 at given surface 
density, and therefore h/r by a factor ~ 1.7, compared to 
solar composition material. 

For completeness, we note that significant outflow may 
accompany accretion flows that are substantially super- 
Eddington. For computational expedience we only consider 
the slim disk in the super-Eddington limit, since the addition 
of outflows would introduce additional free parameters. How- 
ever, given the strong potential for outflow, we stress that our 
accreted masses represent lower limits as to the amount of 
mass that may actually be accreted during segments II and 
in. Similarly, the inferred accreted masses shown in Figure 1 
may significantly underestimate the true remnant masses left 
over after the coUapsar event. 

CG09 present analytical models for the afterglow light 
curves powered by fall-back disks. The main strength of 
model is the universal decay law dlogLacc/rflogt — —1.3 
characterizing fall-back disks without the external, confining 
tidal torque of a companion star (Cannizzo, Lee, & Goodman 
1990). This law seems consistent with the late time decay of 
GRBs seen in X-rays. The early time decay, in particular the 
steep decay for t < 10"^ s, following the prompt emission, 
is difficult to explain within the fall-back disk scenario and 
may be giving us information about the radial density pro- 
file within the progenitor (Kumar et al. 2008b). In this work 
we present results of time dependent numerical calculations 
using a general accretion disk code (Cannizzo et al. 2010), 
which uses the input physics detailed in CG09. The boundary 
conditions are as follows: (i) at the inner disk edge, taken to 
be 10^ cm for sGRBs and 10^ cm for IGRBs, matter is in- 
stantly removed as it arrives, while (ii) the outer grid point is 
placed at such a large radius, typically 10^'' cm, that during 
the course of the run the outer edge of the spreading fall-back 
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Fig. 2. — Model light curve showing the accretion-derived X-ray flux from 
the fall-back disk in a IGRB, taking Mbh = SM© and A/disk = lO""'*/©. 
Also shown is the Swift XRT light curve (Evans et al. 2007, 2009) of GRB 
060729 (in blue). Given the value Airdl = 1.15 X 10^^ cm^ for GRB 
060729 (z = 0.54), the efficiencies assumed in the two scaled model light 
curves (in red) to convert from rate of accretion onto the central engine to 
XRT fluxes are ft~^^ = 8.8 X 10"^ (upper) and 1.4 X 10"^ (lower). The 
accretion rate values given on the right hand side axis correspond to an overall 
efficiency /e"^)- = 4.9 X 10"'^. The five flux values at late times 10^ s 
< t < 10* s (in green) come from Chandra observations (Grupe et al. 2010). 
The last data point lies at i — TO = 642 d. The initial profile is determined 
from Heger progenitor model 12SE as explained in the text. The flat decay 
portion in the model light curve corresponds to the time during which the 
initial 2{r) is being redistributed into an accretion disk with a self-consistent 
radial profile. The lower panel shows the evolution of the locally defined rate 
ofdecayo = — d log Lx/d log t. The two quasi-constant portions, a ~ 0.4 
and 1.2, are indicated by the heavy dashed red lines. The long term decay rate 
~ 1.2 is well-characterized by the analytical solution for electron scattering, 
gas pressure disks, a = 19/16 (Cannizzo, Lee. & Goodman 1990). 

disk never reaches it. The large dynamic range in disk radii 
''outer/nnnor necessitates ^ 1500 — 3000 grid points. In ad- 
dition, no fresh material is added during a run. Thus the evo- 
lution is set entirely by gradients within the disk, unlike the 
standard Shakura-Sunyaev disk fed at a constant rate in the 
outer edge which approaches a steady-state M{r) = constant 
with time. 

4. LORES 

LGRBs are thought to be associated with the explosions 
of massive stars (MacFadyen & Woosley 1999, Woosley & 
Heger 2006). The initial failure of the SN is the reason 
for the GRB in the collapsar model. The hosts for IGRBs 
tend to be subluminous, irregular galaxies rich in star for- 
mation (Fruchter et al. 2006). If the long-term X-ray light 
curves of GRBs are indicative of feeding from a fall-back ac- 
cretion disk, then from Figure 1 we see that no more than 
^ 10^"* — IO^'^A/q of material in the fall-back disk is 
needed, for nominal assumptions about the efficiencies. In 
other words, for a ^ IQMq progenitor, a mass fraction only 
~ 10"^ — 10^"* (excluding the ^ IOA/0 which ends up within 
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Fig. 3. — The disk evolution accompanying Figure 1 . The upper panels 
show S(r) (solid lines, top to bottom) at t = 5, 25, 100, 300, and 600 d. 
The dashed line indicates the initial S(r) distribution, and the parallel dotted 
lines indicate the critical values associated with the dwarf nova limit cycle 
which occurs at the transition between neutral and ionized gas. The lower 
panels show the evolution of h/r corresponding to the same times as in the 
upper panels. 
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Fig. 4. — The evolution of the boundary between the thick and thin disk 
states, described in the text, for the run shown in Figures 2 and 3. 

t < 10 — 10^ s in the BH) is required to survive and persist 
in the vicinity of the progenitor to power the long-term light 
curve. If a significant fraction of the mass that tries to ac- 
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Crete onto the BH is ejected, the remnant mass around the BH 
could be substantially larger. We stress that in this discus- 
sion we only include mass accretion beginning in segment II, 
since segment I is too steep to be accounted for easily in the 
accretion disk scenario. 

Figure 2 shows a light curve for A/bh = 5AIq and A/disk = 
10^'' Mo, overlaid with XRT data from GRB 060729. The 
initial E(r) profile is taken from Heger progenitor^ 12SE 
(Woosley & Heger 2006) in the following way: The progen- 
itor is collapsed in cylindrical coordinates along the z-axis to 
obtain Sinitiai(?'), and then scaled down in mass with a radial 
factor c>c (r/R^) so that J 27rr(irl]iiiitiai(7') equals 10^^ Af©. 
A radial scaling is also applied so that the integrated specific 
angular momentum in the final state, taken to represent the 
disk and therefore to be Keplerian, equals that in the progen- 
itor This effectively introduces a radial scaling ^/^totZ^crit 
to the Einitiai(T') distribution. Due to the very long plateau 
of GRB 060729, we adopted a progenitor model with a large 
radius. (The value logtj*j(s) ~ 4.7 for GRB 060729 is at the 
upper end of the distribution of log t jj values given in Dain- 
otti et al. 2010.) The lower panel of Figure 2 indicates the 
local value of the temporal decay index a associated with 
the model. There are two power law segments evident in the 
model decay light curve, a brief one a ~ 0.4 between 10'' and 
10"' s, and a longer one a ~ 1.2 between 10^ and 10^ s. The 
latter value a ~ 1.2 is close to that expected from analyti- 
cal models, 19/16 (Cannizzo, Lee, Goodman 1990). Figure 3 
gives the evolution of surface density I](r, t) and local aspect 
ratio h/r, and Figure 4 indicates the transition point between 
the slim disk and thin disk states, i.e., h/r dropping below 
(3/4)1/2. 

The steep decay associated with segment I is not accounted 
for in our model, and may well be due to a separate physical 
process, such as the accretion of the progenitor core (Kumar 
et al. 2008b; Lindner et al. 2010, see their Fig. 2). The time 
scale til for the plateau associated with segment II comes out 
naturally in the models, given the ^ 10^^ cm radius for the 
progenitor — assuming that the initial fall-back mass distribu- 
tion roughly traces the envelope of the progenitor (after tak- 
ing into account the radial scaling factor ■\/jtot /jcrit)- The 
observed X-ray decay for GRB 060729, the GRB which has 
the been observed for the longest time in X-rays (Grupe et 
al. 2007, 2010), does not match precisely that of the fall-back 
disk. There may be various systematic effects that could ac- 
count for the difference. The most obvious effect, a variable 
cosmological K-correction concomitant with the fading in X- 
ray flux over a dynamic range of ^ 6 decades, cannot play a 
strong role, as the spectral index does not vary significantly 
over the long term from its mean value ~ 2. There may be 
some time-variable physics associated with the X-ray emis- 
sion from within the jet that could affect either enet or / ver- 
sus time. Our calculations only represent the accretion disk 
which powers the jet and we take enct = / = for simplicity. 
Modeling the jet emission is beyond the scope of this work. 

At late times 10^ s < ^ < 10^ s there is a deviation from 
the canonical power law decay because E(r'outer) drops to the 
point at which the dwarf nova limit cycle becomes active (in- 
dicated by the parallel dashed lines in Figure 3), and a cooling 
front is launched from the outer edge that instigates a transi- 
tion from ionized to neutral gas. This quenches the source of 
accretion onto the central engine. The late time Chandra ob- 



servations of GRB 060729 (Grupe et al. 2010) appear to co- 
incide with the timing of this drop-off, which may therefore 
represent a cooling transition front in the outer disk which di- 
minishes the rate of accretion onto the central engine, rather 
than a jet break. 

5. SGRBS 

SGRBs are thought to be caused by the merger of two neu- 
tron stars (Eichler et al. 1989; Paczyhski 1991; Narayan, Pi- 
ran, & Kumar 2001; Rosswog & Ramirez-Ruiz 2002). The 
hosts for sGRBs are indicative of an older population with 
less active star formation (Fox et al. 2005; Gehrels et al. 2005; 
Villasenor et al. 2005; Bloom et al. 2006), and also tend to 
be ^ 5 times more spatially extended than those of IGRBs, 
matching the ^ 5 times wider projected spatial distribution 
of sGRBs vs. IGRBs within their hosts (Fong, Berger, & Fox 
2010). Due to their faintness relative to IGRBs, sGRBs are 
much less well-studied in X-rays. Indeed, there is only one 
certifiable example of a sGRB which has a long and well- 
sampled X-ray light curve, to the extent that a statement can 
be made as to the presence of a plateau - GRB 051221A. 




r(cm) 



r(cm) 



' models available at http : //2sn . org/GRB2 



Fig. 5. — Model light curve showing the accretion-derived jet power from 
the fall-back disk in a sGRB, taking Mbh = 3Mq andAfdisk = 10~^A/q. 
Also shown is the Swift XRT light curve (Evans et al. 2007, 2009) of GRB 
051221A (in blue). Given the value Andf = 1.18 X 10^^ cm^ for GRB 
051221A (2 = 0.5465), the efficiencies assumed in the two scaled model 
light curves (in red) to convert from rate of accretion onto the central engine 
to XRT fluxes are = 2.4 X 10"^ (upper) and 4.3 X 10"^ (lower). One 

sees fairly rapid decay initially due to the presence of a significant amount of 
accreting material at small radii. The brief shoulder in the light curve is due 
to the delayed accretion of the spiral arm of ejected NS matter at large radii, 
motivated by detailed SPH calculations (Rosswog 2007). The bottom panels 
show the surface density evolution, at equally spaced intervals 5t = 2 d. The 
dashed line indicates the initial E(r) distribution, and the parallel dotted lines 
indicate the critical values associated with the dwarf nova limit cycle. 
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The importance of the outward spreading of the disk formed 
from the NS-NS merger was clearly demonstrated by Met- 
zger et al. (2008a, see their Fig. 1). The initial radial density 
profile, and hence S(r) profile after after angular momentum 
conservation has vertically compacted the merger remnants 
into their accretion plane, is of course much more radially 
condensed than that expected for IGRBs. Theoretical guid- 
ance on a potential S(r) profile to begin the calculations, par- 
ticularly at large radii, is scarce, yet there are indications of 
small amounts of matter ~ 10~^ — 10~^ Mq ejected at early 
times that may have high angular momentum, and therefore 
circularize at large radii (Rosswog 2007). One potentially im- 
portant effect we cannot model is the disruption of a portion 
of the disk by strong nucleosynthesis in the NS supphed ma- 
terial as it expands rapidly to subnuclear densities (Metzger 
et al. 2008a, 2008b, 2010). Lee, Ramirez-Ruiz, & Lopez- 
Camara (2009) find that strong winds can be launched from 
the surface of post-merger NS-NS disks, powered by the re- 
combination of free nucleons into a— particles. 

Figure 5 shows the evolution of a fall-back disk of poten- 
tial relevance for the aftermath of a NS-NS merger. We take 
Mbh = SMq and Mdisk = 10"'"' Mq. Unlike the much 
longer evolution shown in Figure 2 for IGRBs, in this case 
the outer edge of the disk is still freely expanding to larger 
radii by the end of the run (indicated by the narrow spike in 
S(r) between 10^^ cm and 3 x 10^^ cm shown in the second 
lower panel). The brief plateau at ~ 10^ s results from the 
ad hoc introduction of the small amount ^ 10^^ Mq of high 
angular momentum material at the large circularization radius 
10^° cm. The high efficiencies associated with BH accretion 
show the potential for a small amount of material <C IMq to 
have a dramatic effect on the long-term light curves as regards 
brief plateaus or inflection points. 

6. POTENTIAL APPLICATION TO A DAINOTTI-LIKE RELATION 

In Section 2 we show that the Dainotti relation Ljj cx i*n^ 
as originally envisioned may be due to the observational bias 
against detecting and characterizing faint plateaus. Neverthe- 
less, the accreted mass estimates SM inferred using the Ljj 
and fjj values from Dainotti et al. (2010) are interesting in 
the context of this work, and one might legitimately ask what 
theoretical prediction for the Ljj (tjj) relation the fall-back ac- 
cretion hypothesis would make, given a hypothetical physical 
constraint in which 6M were held constant and the initial ra- 
dius of the fall-back disk were aUowed to vary. Figure 6 shows 
the results of seven runs in which we fix SM = 10~^Mq and 
vary the initial outer radius of the fall back disk. One can see 
a clear inverse relation between the duration of the plateau ijj 
and the luminosity at the end of the plateau Ljj. The low— 2; 
behavior of 5M seen in Figure 1, however, indicates that a 
spread in SM may be more realistic, in which case one would 
not expect to be able to use the theoretical prediction of the 
Ljj(tjj) relation as a useful discriminant for the theory. 

7. DISCUSSION AND CONCLUSION 

We have presented time dependent calculations of the fall- 
back disk scenario to account for the long-term X-ray light 
curves for GRBs. For IGRBs, an initial radial scale of a stel- 
lar radius ^ 10^^ cm gives a natural viscous evolution time 
of ^ iC s to redistribute the matter into a quasi-steady disk, 
roughly consistent with the observed plateau duration. GRB 
060729 which we model in the context of a IGRB had one 
of the longest plateaus ever observed, thus more generally 



10" 




Fig. 6. — Model light curves for IGRB parameters, keeping the initial 
fall-back disk mass constant at 10~'^Mq but varying the initial radius and 
overall normalization. For the central light curve ro = 7x 10^° cm. For 
each successive light curve ro is increased a factor of two going to the right, 
and decreased a factor of two going to the left. To keep SM constant, the 
overall normalization on E(0) is varied by an additional factor of 2.5 for 
each successive run. 

we anticipate progenitor radii, or more specifically, circular- 
ization radii, given by ~ i/jtot/jcrit-R*. to He in the range 
~ IQio - 10" cm. The rate of decay for GRB 060729 is 
close but does not match the models precisely, which may 
hint at time- variable emission processes that would affect our 
adopted efficiencies fnci and /. It is interesting to note that 
our effective values of the efficiencies, accretion plus beam- 
ing, required to match the observed X-ray flux levels are com- 
parable between IGRBs and sGRBs: fe~^^ ~ IQ-^ for GRB 



060729 and /e,^ 



3 X 10-2 for GRB 051221 A. If the accre- 



tion efficiency end is about the same for IGRBs and sGRBs, 
this may indicate that sGRBs are less beamed by a factor of 
^ 3 compared to IGRBs, roughly in line with previous re- 
sults (Watson et al. 2006, Grupe et al. 2006). This sinularity 
between IGRBs and sGRBs afterglows is not entirely unex- 
pected, given the general similarities in their afterglow prop- 
erties (Gehrels et al. 2008; Nysewander, Fruchter, & Pe'er 
2009). 

The concept of powering the long term X-ray light curve of 
GRBs by accretion onto the BH represents a departure from 
the standard model in which the fading corresponds to the de- 
celeration of a baryonic jet. In the accretion model, the jet 
itself would be very Ught, perhaps composed almost exclu- 
sively of Poynting flux, and the long term fading would not be 
due to variations in the Lorentz factor and beaming factor. Re- 
cent high fidelity GRMHD calculations of BH accretion sup- 
port the idea of a high Lorentz factor jet with minimal baryon 
loading (McKinney & Narayan 2007ab). The accretion sce- 
nario involving fall-back of the progenitor core has been ex- 
amined in detail by Kumar et al. (2008a, 2008b) and Lindner 
et al. (2010). Kumar et al. (2008b) argue that constraints 
may be placed on the density profiles and radii of the progen- 
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itor core and envelope, as well as their rotation rates relative 
to break-up. Lindner et al. (2010) examined the progenitor 
core fall-back scenario in much greater detail with the adap- 
tive mesh refinement FLASH code, where their calculations 
are done in cylindrical coordinates. Starting with a Heger pro- 
genitor, they follow the fall-back evolution of the progenitor. 
Their Fig. 2 shows the potential for the model to obtain a 
steeply decaying light curve as is seen in X-rays for segment 
I, where a ~ 3 in the Zhang et al. (2006) schematic. In their 
Discussion, Lindner et al. mention several caveats, such as 
the lack of nuclear physics in the inner disk, the neglect of 
the MRI, and the lack of modeling the axial relativistic jet. A 
more basic concern is simply the range of rate of accretion 
in comparison to that inferred from observations. Fig. 2 of 
Lindner et al. shows that if progenitor core fall-back is the 
correct explanation for segment I in the X-ray light curve, the 
rate of accretion onto the BH during this phase varies between 
about 0. 1 and lO^'''Af0 s^^, whereas from Fig. 2 of this work 
we see that, for nominal assumptions about the accretion ef- 
ficiency, the rate inferred from observations on segment I for 
GRB 060729 varies between about 10 and 3 x IQ-^M© yr-\ 
lower by a factor ^ 10'"' — 10® than in Lindner et al. (2010). 
One does in fact expect for the theoretical accretion rate in this 
context to be an overestimate, given the potential for outflow, 
but the discrepancy here seems extreme. There are several 
possibiUties for this discrepancy. The vertical scale shown in 
Fig. 2 of this work is already super-Eddington by between 
about one and eight orders of magnitude, and therefore the 
range in the Lindner et al. calculations would be correspond- 
ingly greater. It may be that even highly advective disks can- 
not accommodate that much accretion onto the BH, and that 
a significant fraction of the material is blown away before it 
can accrete. This would have to occur, however, in such a way 
that it did not interfere with the propagation of the jet. Also, 
for the shape of the decay light curve calculated by Lindner et 
al. to correspond to segment I, the ratio of accreted to ejected 
gas would have to remain constant. If it varied significantly, 
the accretion-derived luminosity would have a different decay 
power law. Another possibility is that the Lindner et al. cal- 
culations greatly over-estimate the fall-back, in which case a 
much greater fraction of the progenitor would be ejected on a 
time scale < 10^ s. Lindner et al. note that they do not find 
evidence for the thin disk hypothesized by CG09 as charac- 
terizing the fate of the progenitor envelope fall-back. Given 
the apparent mismatch in accretion rates between their theory 
and the observations, and the findings of this work, their crit- 
icism could have at least two mitigating factors: (1) In Figure 
4 we see that for the period covered by the Lindner et al. cal- 
culations, namely t < 10'^ s, a significant fraction of the disk 
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lies within rtians, i-c, the more spatially extended slim disk 
rather than the thin disk. (2) Given the apparent mismatch 
between theory and observation for the rates of accretion, the 
actual densities within the volume formerly occupied by the 
progenitor may be much less than calculated by Lindner et al. 
(2010), which would interfere far less with an accretion disk 
formed from fall-back debris. 

At late times Chandra observations indicate a steepening 
in the rate of decay (Grupe et al. 2010), which appears to be 
consistent with the onset of a cooling front in the disk. This 
would be an alternative to the standard jet-break interpretation 
discussed by Grupe et al. (2010). For sGRBs, the picture is 
less clear, given that we only have a single well-studied ex- 
ample. We have shown that the presence of even a very smaU 
amount of high angular momentum gas ~ 1O~^M0 can give 
a slight inflection to the X-ray decay, as was observed in GRB 
051221A. As for the overall X-ray light curve, segments 11 
and 111, if only ^ 10^"'' — 1O^^M0 of gas survives either the 
hypernova (IGRBs) or NS-NS merger (sGRBs), then the ac- 
cretion resulting from the ensuing fall-back disk should power 
a long-term jet. 

Lastly, we have shown that the Dainotti relation LJj oc <*n^ 
may be due to an observational bias against detecting and 
characterizing faint plateaus: the relation is governed by 
GRBs at 2; > 1 .5 for which we only detect the upper envelope 
of a broad distribution. Nevertheless, the existence of an ap- 
parent upper limit to the total X-ray energies inferred from the 
X-ray fluences, and therefore the accreted masses if one as- 
sumes accretion onto the central engine as the long term pow- 
erhouse for the X-ray flux, is extremely interesting. For nom- 
inal values of the accretion efficiency enot and the beaming 
factor /, we find an upper limit ~ 10~^ — 10^^ Mq for the ac- 
creted mass during segments 11 and 111. (The lower end of the 
distribution of accreted masses, which is partially revealed for 
GRBs atz < 1.5, may extend down to ~ 10~^ - 10"'^ M©.) 
This means that for a progenitor mass ^ IOM0, only a maxi- 
mum mass fraction ^ 10^"'' — 10^^ of the progenitor survives 
in the vicinity of the progenitor to be accreted as a fall-back 
disk (excluding the ~ lOM© that ends up in the BH during 
the prompt emission and subsequent segment I consisting of 
the steep-decay). This has important ramifications for the en- 
ergetics associated with the hypernova explosion and subse- 
quent removal of most of the progenitor envelope. 

We acknowledge useful conversations with Maria Dainotti, 
Dirk Grupe, Stephan Rosswog, and Brad Schaefer. This work 
made use of data supplied by the UK Swift Science Data Cen- 
tre at the University of Leicester. 
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